The tight-coupling approximation for baryon acoustic oscillations 
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The tight-coupling approximation (TCA) used to describe the early dynamics of the baryons- 
photons system is systematically built to higher orders in the inverse of the interaction rate. This 
expansion can be either used to grasp the physical effects by deriving simple analytic solutions or 
to obtain a form of the system which is stable numerically at early times. In linear cosmological 
perturbations, we estimate numerically its precision, and we discuss the implications for the baryons 
acoustic oscillations. The TCA can be extended to the second order cosmological perturbations, 
and in particular we recover that vorticity is not generated at lowest order of this expansion. 

PACS numbers: 98.80.-k 



The baryon acoustic oscillations (BAO) generated 
when baryons and photons were highly coupled have now 
been observed in the cosmic microwave background [l], Q 
and in the large scale structures [jj. The shape and am- 
plitude of these oscillations cannot be obtained analyti- 
cally and one resorts either to a numerical resolution or 
to a less accurate WKB approximation on the full sys- 
tem of dynamical equations for the cosmic fluids and the 
metric. On the one hand the numerical integration has 
stability issues due to strong restoring forces, and on the 
other hand the precision of the WKB approximation is 
limited. The tight-coupling approximation, which is an 
expansion in the interaction rate [see Eq. (|17[) below], 
is a way to avoid these two issues 0, The orders 
of the approximation are denoted by TCA-n, and one 
can always consider the equations up to TCA-n with n 
sufficiently large to ensure the required accuracy, in the 
range of time for which the expansion is converging. In 
this paper, we first recast the baryons-photons dynamics 
into a total fluid system, and we then present the TCA 
expansion. We also estimate numerically the convergence 
before recombination. A closed form at TCA-1 for pho- 
tons density perturbations is then derived from which 
the main features of the dynamics are deduced. At sec- 
ond order in cosmological perturbations, the TCA-0 is 
sufficient at early times given the precision required for 
all practical purposes. An exception arises for vorticity 
since it is not generated for a single perfect fluid, and we 
recover the known result that it is not generated below 
TCA-1 . Our analysis is of interest for high precision 
BAO computations but also sheds some light on the fluid 
approximation of the baryons-photons system in general 
relativity and cosmology by offering another method for 
solving cosmological perturbations. 



Perturbations 

Geometry: We perturb the metric in the conformal 
Newtonian gauge according to 

ds 2 = a 2 [(1 - 2^)5 13 dx l dx 3 - (1 + 2$)d? ? 2 + 25 i dx i dn] 



where a is the scale factor and r\ the conformal time. Si 
is a vector type peturbation (diS 1 = 0) considered only 
at second order in perturbations since vector perturba- 
tions decay at first order if not sourced. We also do not 
consider the tensor perturbations. A prime denotes a 
derivative with respect to r\, and fi = a' /a . 

Baryons-photons system: The stress energy tensor of 
a species labelled by s is decomposed as 



rr = (ps + P s )« + P s g^ + P s nr 



(1) 



The energy density is perturbed according to 
Ps = Ps(l + $s), velocities are perturbed according 
to u l s = l/a[(l + ^)vl — S 1 ] such that vl matches the 
tetrad components of the velocity We assume that 
baryons have no anisotropic stress (n^ = 0), have 
equation of state Wb = Pb/ph = 0, whereas photons have 
equation of state w 1 = 1/3. The adiabatic speed of 
sound c 2 = P s '/Ps is constant for these two fluids and sat- 
= i« 7 /b- In general for a perfect fluid, if w s is 



isnes c 2 /b 



not constant, c 2 ^ w s and w' s = —3H(c 2 — w s )(l + w s ). 
For photons, free-streaming produces an anisotropic 
stress H]j and in general we have to consider the mo- 
ments of the temperature O^...^, and also the electric 
and magnetic type multipoles (-E 1 ^...^ and i^...^) to 
describe linear polarization. The lowest multipoles are 
related to the fluid quantities up to second order by || 



e, 



Pi 



4p 7 



<« 3) 



(2) 



At first order, there arc only scalar perturba- 
tions, the magnetic type multipoles vanish and 
the temperature scalar multipoles are defined by 
e n ... le =(~l) e (2e-l)l/e\d {ll ...d le) Q l (and similarly 
for electric multipoles) where (...) means symmetric 
traceless part of the indices. In Fourier space, for a given 
mode k with k = |k|, we define dimensionless multipoles 
from the previous ones by O^ = k Q e , and similarly 
for electric multipoles. For velocities, we use instead 
u| = diV s and v s = v s /k, that is w 7 = — @±. 
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Baryons-photons system and tight-coupling 

Fluid equations: In Fourier space, the conservation and 
the Euler equations at linear order are, with s = b, 7 



a = 



1 + w s 



fcu s - 3*' = 



(3) 



kr 2 Ok 

E s = v' s + (1 - 3c 2 s )Hv s + M> + —^S s - =-O s 2 = C s . 

l + w s 5 

The collision terms for baryons and photons are obtained 
from the Boltzmann equation [f| and are given by 



C 7 = r'(ub — Uy) , Cb = — C 7 /i? 



(4) 



with R = 3pb/(4p 7 ), and the interaction rate is 
t' = an e o~T where n e is the number density of free elec- 
trons and or the Thomson cross-section. 

The plasma of photons and baryons is an unperfect 
fluid (labelled by pi) whose energy density and velocity 
are given by 

Ppi = ft . + ^pO<i = E ^ + p *k- ( 5 ) 

s— b,7 s— b,7 

We can infer easily the equation of state of the plasma, 
and by considering its time evolution we find its adiabatic 
speed of sound. They read 



w p i = (3 + 4R) -1 , c 2 pl = [3(l + R)Y 



(6) 



If we define a reduced energy density contrast by 
(1 + w s )A s = (1 + w p i)S s , we obtain from Eqs. (0 the 
very simple first order relations 



S p i = Bar(A b , A 7 ), v p \ = Bar(u b , v 7 ) 



(7) 



where we defined Bar(A b , A 7 ) = (RX b + X y )/(1 + R). 
The conservation equations are then rewritten in a com- 
pact form as 



a 



A. 



1 + Wpl 



- kv^ - 3*' = . 



(8) 



In order to fully characterize this two fluids system we 
also need to describe the differences between them by the 
entropy perturbation S = 3[(pb/pb) 1 ^ 3 — {p-y/P-y) 1 L 
and the velocity slip = u£ — n't. Both quantities van- 
ish for adiabatic initial conditions at any order in pertur- 
bations. We obtain at linear order 



(l + u;pl)S = A b 



V = «b — v~ 



(9) 



In the rest of this paper, we abbreviate w p \ to w and c 2 \ 
to c 2 . The photons-baryons system is either described 
with the intrinsic fluid variables (<5 7 , <5b, i>b, w 7 ) or with 
the total fluid variables (S p i, v p i, S, V) but the latter 
choice of variables is better suited for the TCA. 



Tight-coupling expansion: From the definition ([T]). 
considering Bar(Cb,C 7 ) and Bar(Eb,E 7 ) with Eqs. (|5J) 
and (@| leads to the plasma equations 



C p i = U 
E p i = k 



RS 



(l + R) 2 
RS 



20i 



3(1 + R) 2 {l + R) 5 



(10) 
(11) 



Taking the differences E b — E 7 and Cb — C 7 , we also 
obtain the entropy and velocity slip equations 

S'-kV = 0, (12) 
-(1 + R)/Rt'V = 

V + n{v p i + 3c 2 V) + k(^Q 2 



(13) 



3(1 + w) 



Rc z S 



The dynamics of the quadrupole O2 is inferred from the 
Boltzmann hierarchy i|. Combined with the dynamics 
of Ei we get 



r'9 2 = 2E', 



8v 7 
~9~ 



49 3 



10g 3 
21 



(14) 



where we must then use that i> 7 = v p \ — RV/(1 + R). 
The same method leads to 



t'E 2 



-3EL 



03 

7 



5E 3 
7 



(15) 



For higher order moments, the Boltzmann hierarchy 
leads directly to 



r'E e = -E't + k 



(I -!)(* + 3) 
{£ + l){2£ + 3) 



Ei+i — 



21- 1 



-E. 



21- 1 



-0/ 



(16) 



When combined with the Einstein equations to determine 
the perturbations of the metric, these equations com- 
pletely determine the system. However, at early times 
t' oc a~ 2 and the restoring forces can be huge and require 
very small numerical steps. The TCA consists in not 
solving the dynamical equations for (V ,0£>2,-E£>2) but 
to use instead the expressions (|13M16p multiplied by 1/t' 
to obtain expressions for these variables as constraints 
which arc functions of plasma variables but also func- 
tions of themselves. To obtain a closed system out of 
this infinite recursion, all variables are expanded in pow- 
ers of the TCA parameter e in the form 



X 



p=0 



with 



n/r' , (17) 



and we use that F (0) = and G$ 2 = E}% = for 
n < (I — 2). The expansion (|TT|) can be obtained by 
replacing recursively the constraints in themselves, and 
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getting rid of the time derivative by using the plasma 
and entropy equations. We stop at the desired TCA or- 
der, gaining one order at each recursion, and then the 
constraints are replaced in the plasma and entropy equa- 
tions. We obtain finally a closed and first order differen- 
tial system up to the desired TCA order. More precisely, 
the result is stated as 

(*): Defining the dynamical and constrained variables 
Y = (v ph 6 ph S,$,V) and Z=(V,e 2 ,E 2 ,e 3 ,E 3 ,...), 
then to solve for their dynamics there exists two sets of 
matrices M p and N p , the coefficients of which are sums of 
products of the (k/H) q , the n" {q+1) d q H /dni } and some 
functions of R (this type of matrix is called generic in 
the following), such that the equations to solve are 



p=0 



Z = ^VN p .Y (18) 

p=0 



with No = 0. This implies in particular that there exists 
generic matrices K p such that Z' = H Yl^Li e p K p .y 
[obtained from K p ee (e p N p y/{He p ) + g = i N r - M P -r]- 
If we define the variables up to TCA-n by 

X <n ee YJp=o eVx(v \ thcn Y<n and Z <n are obtained 
from a truncation of the expansion (|18|) . 

In order to show (*) we will determine these matri- 
ces recursively. Hereafter P, Q, U, W, B, and C are 
generic matrices. First, the dynamics of \P from Einstein 
equations is of the form iff' = TLP.Y, but <!> is found 
from a constraint $ — cx 02 and its dynamics depends 
on 6j [H. Using this, the plasma-entropy equations are 
of the form Y' = H[Q.Y + V.Z + W.Z'/H], the depen- 
dence in Z' being just on Oj'. The M p are not indepen- 
dent since M = Q and M p = U.N p + W.K p if p > 0. 
Hence, the K p and M p are deduced from the N p . For 
p = 0, M = Q and N = K = 0. At TCA-0, it is as if 
there was a single perfect fluid. Now, if the matrices are 
known up to a given p > 0, that is if we know Y <p and 
Z <p (and thus Z <p ), then we can find K p+ i, and using 
that the constraints for Z is of the form 



Z 



A.Z'/H + B.Z + C.Y 



(19) 



where A is a constant matrix, we replace Y <p , Z <p and 
Z <p in the right hand side. We then get N p+ i from 



N 



P+i 



A.K„ + B.N, 



(20) 



for p > 1 and Ni = C otherwise, from which we find 
M p+ i. At early times (r')' = -2Ut' , and R' = UR, and 
this is why matrices are of generic type. For a realistic 
case, neutrinos an cold dark matter should be considered 
and added to Y . Note also that when recombination 
starts, t' is no more scaling like oc a~ 2 and one must use 
its complete expression. However, (r')' factors appear 
only at TCA-2. The public code CAMB uses the 
TCA-1 with intrinsic fluid variables instead of total fluid 



variables, and such factors appear already at that order. 
With this algorithm, we can find the equations up to a 
given TCA-n if one uses abstract calculus to perform the 
recursion. Our approach clarifies the method adopted 
in Ref. [lj| where (*) is implicitely used up to TCA-2. 
At early times, it can be used to avoid instablities due to 
the high interaction rate. Using higher orders of the TCA 
enables to improve the accuracy, given that up to TCA- 
n th precision is roughly of order (k/He) n = (k/r') 71 for 
large modes and e n otherwise. In Fig. [T] we plot the error 
between the full numerical integration and the successive 
orders of the TCA for v 1 . 

v y and |v r TCA -i> r | TCA=0,1,2 




-1.0 -0.5 
Log l0 [4/3 Rj 

FIG. 1: Top thick line: exact v 1 with k = 0.2Mpc~\ The 
differences with the solutions up to TCA-0, 1, 2 are in continu- 
ous, dashed and dotted lines. The orders expected (k/r') 1 ' 2 ' 3 , 
are depicted in thin lines of the same type. 

Baryon acoustic oscillations: Qualitative features of 
the BAO can be inferred by computing the evolution of 
i5 7 at TCA-1. Considering the combination C 7 + fcE 7 
and expressing 02 up to TCA-1 and V up to TCA-2 in 
order to obtain V/e up to TCA-1, and further ignoring 
the variations of the potentials, we obtain 



1 



eR 2 



-Si 



(1 
UR 
l + R 



Rf 



k 2 c 2 6. 



16 
15 



1 



R 2 



R 



eR 1 



l + R V 1 + -R 



1 + -R 



= - A -k 2 *. 
3 



where r = (lnr')'/?^ and the dots in the friction term 
represent terms which are suppressed by a factor (R/k) 2 
with respect to the dominant TCA-1 term. We thus re- 
cover that up to TCA-1, the viscosity is damping the 
oscillations [llj, but we also find that the pseudo-period 
of oscillations, and thus the sound horizon, is slightly 
modified by the TCA-1 terms. However for large modes, 
k/r' > 1 before e > 1 and the effect of viscosity is ex- 
pected to be dominant over the modified pseudo-period 
variation. This approach is of limited precision, since the 
TCA cannot be trusted around recombination, and only 
the full numerical integration makes sense. 
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Second order cosmological perturbations: vorticity 

At second order, the same procedure can be followed. 
Given the precision that we need to reach, the TCA-0 is 
sufficient, and the equations are those of a perfect fluid. 
They can be found e.g. in Ref. [IH (with a different 
definition of the velocity perturbation). There is however 
a case where the TCA-0 is not sufficient. Indeed, vorticity 
is not generated for a perfect fluid [l3| and thus not at 
TCA-0 p. One needs then to consider the equations up 
to TCA-1 at least in order to describe the generation of 
vorticity. The Euler equation for a perfect fluid reads 



f 



-diS s + - 4*'v| 



+ (f - 3c 2 s )nv! + 

2 [{S B vty + H(l - 3w 8 )5 B vl + 6 s di$] 
)-(* + #) [v] 
h (1 - 3c 2 s )Hv! 

k 2 ) 



l + w s 

dM4 



H{1 - 

s 

W s 



1 



-3cJK]-ft(* a ) 
d t 5 S +di$ 



(c 2 )' 



6H(l + w s ) : 



(21) 



The vorticity is defined by = _L s jU a ± s /V[ Q ,'U^i where 



\ (Xftv — X 



V fj), and the projector is given by 
There is no vorticity a first order 
since we have discarded the vector modes, but at second 



X \M - 2' 
_L s ,r = 6" + ulu" 



order its expression is 



CLUJ,; 



(22) 



For a collisionless perfect fluid, that is if C^ s = 0, then 
from the Euler equation, we can show that the evolution 
of vorticity is dictated by 



(2 - 3<S)«w? 







(23) 



which implies that [p s (l + w s )a 5 wf J ] is conserved. At 
second order, the velocity of the plasma is given by 

Bar(l + 5 b ,l + * 7 )u^ = Bar[(l + i b X > (l + 5 7 )u^]. (24) 

Using this, we find as expected that at TCA-0 the 
baryons-photons plasma behaves like a perfect fluid, since 
Bar(E^ , Ej' ) = E? 1 , and the plasma vorticity satisfies 



^ (0) = 0. 



■«5 



However in full generality 
R 



[1 + RY 



3(l + w p i) 



o 



with 6™ m = S p i - 3H(1 



Wp\)v p \, and where fijj = if 
for photons we neglect the anisotropic stress and use a 
perfect fluid description. This equation without the last 
term (which is related to the quadrupole generated by the 
mixing of the fluids) matches the expression of Ref. [3] , 
given that the non-adiabatic pressure perturbation is 



-Ppl 



R(l + w p i) 
3(1 + R) 2 



S. (25) 



However we do need to consider the anisotropic stress at 
TCA-(n > 1). For completeness, its contribution is 



1 + c 2 v?(c 2 d k Tl\ ) 

1 + w Jl & 



Rd\ Sd k lV 

3 (1 + J?) f + + 6 ^ + 3 ^ U M 

In this paper, we have formulated the TCA in the total 
fluid variables since it can then be extended easily up to 
any order through a recursions on the equations. Higher 
orders in the TCA can be used to speed up and refine lin- 
ear Boltzmann codes 0, [loj]. At second order, a TCA-0 
solution is generally sufficient, e.g. for the computation 
of non-Gaussianity generated by non-linear effects 15 1. 
However, the vorticity in the baryons-photons fluid is 
generated at least at TCA-1, not only from gradients of 
non-adiabatic pressure perturbations (see Refs. [tl IT3|). 
but also from gradients of the anisotropic stress of pho- 
tons. This becomes relevant for the numerical estimation 
of the seed magnetic field created by vortical currents [l6| 
since it is related to the existence of vorticity. 
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